Local elastic response measured near the colloidal glass transition 
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We examine the response of a dense colloidal suspension to a local force applied by a small magnetic bead. 
For small forces, we find a linear relationship between the force and the displacement, suggesting the medium 
is elastic, even though our colloidal samples macroscopically behave as fluids. We interpret this as a measure 
of the strength of colloidal caging, reflecting the proximity of the samples' volume fractions to the colloidal 
glass transition. The strain field of the colloidal particles surrounding the magnetic probe appears similar 
to that of an isotropic homogeneous elastic medium. When the applied force is removed, the strain relaxes 
as a stretched exponential in time. We introduce a model that suggests this behavior is due to the diffusive 
relaxation of strain in the colloidal sample. 



I. INTRODUCTION 

Glass is an amorphous solid: despite the lack of long- 
range order, a glassy material is elastic rather than vis- 
cous. As a glass- forming material is cooled, its viscosity 
rises dramatically by many orders of magnitude; flow be- 
comes difficult and slow. The elastic behavior of a glass, 
then, could be reframed as the material being probed 
on time scales too quickly for liquid-like flow to occur. 
The origins of this elasticity and the nature of the glass 
transition are a quite active area of researcftt"— . 

A related question is to what extent the macroscopic 
elastic behavior extends to the microscopic scale. For ex- 
ample, simulations have seen that the elastic moduli are 
spatially heterogeneous in some cases^ii. Certainly the 
macroscopic elastic response treats the material as a con- 
tinuum, whereas on a scale of the constituent molecules 
this could be a poor approximation. 

Colloids are a simple model system which can be 
used to study the glass transition and the properties 
of glassy materials^"—. Colloidal suspensions are com- 
posed of solid particles in a liquid. The particles dif- 
fuse due to Brownian motion, but this diffusion is im- 
peded at high particle concentration. In many experi- 
ments, colloidal particles only have short-range repulsive 
interactions, and the particles can be approximated as 
hard sphere o^^'^^ . In such samples the control param- 
eter is the volume fraction (p and glasses are found for 
4> > (j>g O.SSiSiii. Near the glass transition, the viscos- 
ity rises quite dramatically, although some evidence sug- 
gests that perhaps it truly diverges at a point above (j)^- 
In particular, an alternative divergence point is "random 
close packing" (rep). This is the largest volume fraction 
possible for a sample that is still amorphously packed 
(rather than crystalline)-"*. The volume fraction (j)rcp is 
known through simulations and taken to be ~ 0.64, which 
agrees with early experiments done with ball bearings^^. 

Of course, the presence of the continuous liquid sur- 
rounding the colloidal particles is important for under- 
standing the flow properties of colloids. Approaching 
the colloidal glass transition from (p < (pg, samples are 
viscoelastiC"^. Their properties are described by both vis- 



cous and elastic moduli which are frequency-dependent. 
The viscosity mentioned above is understood to be the 
low-frequency limit. Indeed, if a glass is defined as an 
amorphous solid - a sample that does not flow - then 
it is important to recognize that whether or not it flows 
depends on the time scale of observatior^i^. 

While macroscopically one considers viscous and elas- 
tic behavior, microscopically one considers diffusion. A 
molecule or tracer particle in a fluid sample has a nonzero 
diffusion constant, which decreases to zero as the glass 
transition is approached. The decrease of the diffusion 
coefficient is attributed to caging. On short time scales, 
particles diffuse in a "cage" formed by their neighbors. 
On longer time scales, these cages rearrange and particles 
can move throughout the sample. As the glass transition 
is approached, the cage rearrangements occur less fre- 
quently, thus decreasing the diffusivitjJ^"— . These cages 
provide a sort of elasticity for individual particlea^S^^. 
Particles which try to move away from the centers of 
their cage experience a restoring forcei^"— . If a constant 
external force is exerted on a particle, it will slowly move 
through the colloidal suspension as cages rearrange, al- 
though the force needs to be kept small when <p is close 
to (pg to avoid nonlinear behavior—"—. If a sufficiently 
high external force is applied to a particle, it can break 
the cages and move through the sample more freely (typ- 
ically with a nonlinear relationship between the force and 
resulting velocity) ^^i^^'^^~ '— , disturbing and rearranging 
particles as it moves^SiH. 

In this paper we probe the elastic response of a dense 
colloidal suspension by locally exerting a small force on a 
magnetic bead. Our studies are conducted on a time scale 
faster than the sample can relax due to diffusion. We find 
that the sample responds elastically, with a Young's mod- 
ulus E that rises as the glass transition is approached, 
and a Poisson ratio a equal to 1/2. We also study the re- 
laxation of the magnetic bead after it has been displaced 
and the force removed. This relaxation is faster for higher 
volume fraction samples. We present a model that cap- 
tures the stretched exponential character of the relax- 
ation, by assuming that in our colloidal sample stress 
diffuses away to infinity at long times. 
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II. EXPERIMENTAL METHODS 

The colloidal suspensions are made of poly- 
(methylmethacrylate) particles, sterically stabilized 
by a thin layer of poly-12-hydroxystearic acid^. The 
particles have a radius a = 1.55 /im, a polydispersity 
of ^^5%, and are dyed with rhodamine^. The uncer- 
tainty of the mean particle radius is ±0.01 fim.. The 
particles are slightly charged, but their glass transition 
is still at (t)g = 0.58 ± 0.01, signaled by the diffusion 
of particles going to zero on experimental time scales 
(several hours). The colloidal particles are suspended in 
a mixture of cyclohexylbromide/ c is- and trans- decalin 
which nearly matches both the density and the index 
of refraction of the colloidal particles^. The density of 
this solvent is 1.232 g/cw?, the index of refraction is 
1.495, and the viscosity \s rj = 2.18 mPa-s. The particles 
are fluorescently labeled for visualization^. Before 
beginning experiments, we stir the sample with an air 
bubble to break up any pre-existing crystalline regions, 
then wait 20 minutes before taking data. 

A small quantity of superparamagnetic beads (Dynal 
M450, coated with glycidyl ether reactive groups) with 
a radius of 2.25 /xm are added to the colloidal suspen- 
sion. We do not observe attraction or repulsion between 
the colloidal particles and the magnetic beads, in either 
dilute or concentrated samples. The beads are not com- 
pletely monodisperse in their magnetic properties, and 
our calibration finds the variability in the effective mag- 
netic force applied to different beads to be less than 10%. 
Also, the magnetic beads are not density matched, and 
their effective weight is 0.1 pN. This is a factor of ten 
smaller than the smallest horizontally applied magnetic 
forces in our experiments. We study isolated magnetic 
beads at least 35 /xm from the sample chamber boundary 
and from other magnetic beads. 

For some of the data (Figs. I1I2I in particular), we use 
a conventional Leica DMIRB inverted microscope with 
a CCD camera. A magnetic bead appears as a large 
dark circle in our images and its position as a func- 
tion of time is found using standard particle tracking 
techniques^. For these experiments, the typical image 
size is 15 x 40 /im^, with images taken at a rate of 30 
frames per second. The magnetic bead position is re- 
solved to within 0.04 /im; see Fig. [T] for example. 

For other data, we use a ThermoNoran Oz confocal 
microscope. With this microscope, we acquire images 
of area 80 /im x 75 /im at a rate of 30 frames per sec- 
ond (256 X 240 pixels^). The magnetic beads are not 
fluorescent and thus appear black on the background of 
dyed colloidal particles (see Fig. S]). Again, we use parti- 
cle tracking procedures to follow the motion of the mag- 
netic particle with a resolution of 0.04 /im. With either 
the confocal microscope or the video microscope, we only 
track the motion of particles in 2D: for the experiments in 
this paper, the magnetic particle position always remains 
with the imaging plane of the microscope. 

We additionally use the confocal microscope to take 



three-dimensional images of these samples to determine 
the volume fraction'^^ . We count the particles within the 
imaged volume, and convert from the measured number 
density n to volume fraction </) using (j) = (47r/3)na^. 
Due to our uncertainty of the mean value of the particle 
radius a, we have a systematic uncertainty of 2% of our 
(j> values- ''. That is, our (f) values are fairly accurate when 
compared with each other, but a reported value of = 
0.50 is uncertain by ±0.01. 

We use a strong Neodymium permanent magnet 
mounted on a micrometer positioner to control the force 
applied to the magnetic beads. The micrometer accu- 
rately reproduces the magnet position and thus our un- 
certainty in the force between different experiments is 
limited by the magnetic bead variability, rather than the 
magnet positioning. For many experiments, a computer- 
controlled stepper motor attached to the micrometer al- 
lows us to slowly and controUably vary the applied force 
over two orders of magnitude. 

For other experiments, we want to apply a short- 
duration magnetic force. We mount the magnet on a 
linear actuator (Uhra Motion D-A.25-HT17). The mag- 
net is then brought close to the magnetic bead resulting 
in a high force acting on the magnetic bead for a short 
time; the details are discussed below. 

Our experiments are controUed-force experiments, in 
contrast to controlled displacemen t For example, 
this means that particles do not necessarily need to move 
in the direction of the applied force, although we dis- 
cuss below that they appear to always do so. Prior 
work by other groups used laser tweezers to move probe 
particles through colloidal suspensions at a controllable 
velocity^*^^^"^ , finding many interesting results such as 
anisotropy of the disturbed region around the moving 
particle*^ and a decoupling of structural and hydrody- 
namic influences on the particle motion^. These prior 
experiments studied probe particles as they moved over 
long distances, in contrast to our experiments described 
below where the magnetic bead always remains close to 
its equilibrium position. 



III. LINEAR ELASTIC RESPONSE 

In our earlier work, we applied a constant force and ob- 
served the steady-state motion of the magnetic beac]3^. 
For large forces, the velocity of the magnetic bead grew 
nonlinear ly with increasing force, consistent with shear- 
thinning. We also observed that the velocity was es- 
sentially zero below a threshold force. The thresh- 
old force grew dramatically as the glass transition was 
approached^^. 

To test the behavior of the sample below the thresh- 
old force for motion, we vary the magnetic force within 
a range of forces that are below the threshold force for 
motion. To ensure that the forces used are below the 
threshold force for motion, each increase in the force was 
followed by a waiting period of 160 s to observe the sub- 
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FIG. 1. The solid line indicates the applied force as a function 
of time. The points show the measured displacement Ax of 
the magnetic bead. The volume fraction is </> = 0.55. 

sequent motion of the magnetic bead. The applied force 
as a function of time is shown in Fig. [T] by the solid line. 
In this situation, the magnetic bead has a finite displace- 
ment, rather than a finite velocity; the displacement is 
shown by the dots in Fig. [TJ As the figure shows, by ap- 
propriately rescaling the units, the displacement of the 
magnetic bead from its original position is linearly pro- 
portional to the applied force. Additionally, during each 
of the pauses at constant force, the magnetic bead ex- 
hibits slight fluctuations around its equilibrium position. 
Some of this is due to Brownian motion, and some of this 
is due to the uncertainty in identifying the position of the 
magnetic bead (~ 0.04 /im). 

Further evidence for the linearity is seen in Fig. [21 
showing the displacement data plotted as a function of 
the applied force. The data are linearly related, and a fit 
to the data leads to an effective spring constant k — 6.8 
pN//im. The arrows shown in the figure indicate loca- 
tions where the magnetic force was constant. This spring 
constant is quite large, equivalent to 4000 ksT/a^ using 
the radius of the colloidal particles or 8400 kBT/a\i^ 
using the radius of the magnetic bead. 
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FIG. 2. The measured displacement Aa; of the magnetic bead 
as the force is varied; the data correspond to Fig. [T] The 
dashed line is a fit to the data, with the slope leading to an 
effective spring constant k — 6.8±0.1 pN//im, offset vertically 
for clarity. The arrows indicate locations where the force was 
held constant for 160 s (see Fig. [l]). 
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FIG. 3. Applied force as a function of time, for the three 
largest maximum forces (see Table I). For all curves, F{t) — 
for t < s. 



IV. PARTICLE DISPLACEMENT FIELDS 

Unfortunately, the results of Figs. [1] [2] are atypical 
in one important respect: for many experiments last- 
ing longer than 0(100 s), the magnetic bead experiences 
a cage rearrangement and does not return to its origi- 
nal position. The spring constant before and after any 
such displacement is always the same, to within our un- 
certainty. Accordingly, to complement the slow exper- 
iments of Figs. [T] and [21 we conduct experiments with 
intermittent and short pulses of force to see the instanta- 
neous response of the sample on a time scale quicker than 
cage rearrangements. A nondimensional way to consider 
this is the modified Peclet number—. Pe* is the ratio 



of the time it would take for a colloidal particle to dif- 
fuse its own size to the time scale for the perturbation. 
The diffusive time scale is ^400-5000 s for the samples 
we study (0 > 0.4), and the perturbation time scale is 
0.25 s. Thus, Pe* « 2000 - 20000, signifying that Brow- 
nian motion is unimportant on the time scales we study: 
particles do not substantially rearrange their positions 
during our experiment. The recovery to the perturba- 
tion (discussed in Sec.[V| takes 0(10 s), which is still in 
the high Pe limit. 

To discover the origin of the linear restoring force, we 
examine the response of the colloidal particles surround- 
ing the magnetic bead. To produce reproducible initial 
strains, we attach the external permanent magnet to a 
linear actuator as described in Sec. [Til ^-nd move the mag- 
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TABLE I. The five difi'erent maximum forces applied, and the 
integrated impulse I = J F{t)dt. The calibration procedure 
(described in the text) has an intrinsic Fmax uncertainty of 
±0.05 nN and an I uncertainty of ±0.005 nN-s. Due to vari- 
ability between different magnetic beads, for a given magnetic 
bead there is also an overall systematic uncertainty of ±10%. 
Graphs of F{t) for the three largest values of -Fmax are shown 
in Fig. H 
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0.068 


0.75 


0.17 



net toward the sample and then away at maximum speed. 
The force applied is ramped up to a maximum value and 
then just as rapidly reduced. To calibrate the force as 
a function of time, F(t), we use this procedure to exert 
a force on magnetic beads suspended in glycerol. Such 
beads move with velocity v{t), from which we deduce 
the force F{t) using Stokes' Law, F = Q-K-qausv , with 
the viscosity of glycerol -q = 0.934 Pa-s. The resulting 
F{t) data are plotted in Fig.[3]for the three largest Fmax- 
These correspond to the cases where the external magnet 
is moved the closest to the sample. For example, the top 
curve in Fig. [3] takes longer to reach its peak and longer 
to return to = 0, as the magnet has farther to move. 

The peaks of the curves, Fmax, are confirmed by mea- 
suring the velocity of magnetic beads in glycerol while 
the external magnet is fixed in its closest distance to the 
sample for a given forcing protocol. Those results agree 
quite well with the values measured from the F{t) data. 
In the results that follow, we refer to the different F{t) by 
their maximum values -Fmax which are listed in Table HI 
An additional way to quantify the F{t) is by their time 
integral, / — J F(t)dt, yielding an impulse which is ap- 
plied to the magnetic bead. These values are also listed 
in Table m In each case, the ratio //F„iax — 0.25 ± 0.03 
s, suggesting that our choice of using F^ax is correctly 
representing I as well, and that the effective pulse du- 
ration is a quarter of a second. This time scale is short 
compared to the Brownian time scale a^g/DMB = HO s. 

The response of the sample to a pulse is shown in 
Fig.m^b). This is a difference image formed by subtract- 
ing the raw image before the pulse [such as Fig.HJ^a)] from 
the raw image after the pulse. In this case, the magnetic 
particle has moved to the left, as indicated by the white 
crescent on its left side. Because the magnetic bead is 
black and the colloids are white, the colloidal motion is 
indicated by the direction of the black crescents, and is 
also clearly leftward. Furthermore, the overall disturbed 
region of colloids is a fairly smooth function of space. 
Adjacent colloids move similar distances in Fig. Sljb). 

This displacement field is highly reproducible, as is 
shown by creating a difference image in Fig. \^c) between 
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FIG. 4. (a) Raw image of particles, before the force is ap- 
plied, (b) Difference between "before" and "after" a single 
force pulse is applied, (c) Difference between two "after" im- 
ages for two subsequent pulses, (d) As the images are repro- 
ducible, a sequence of eight "before" pictures are averaged 
together, and likewise eight "after" pictures. This picture 
is the difference between these average images. For all pic- 
tures, the scale bar is 10 microns long, the volume fraction is 
4> = 0.49, and the applied force is -Fmax = 0.29 nN. 



two images both taken when the colloids are maximally 
displaced from their equilibrium positions. The differ- 
ence image is nearly completely gray, showing that the 
displacements are virtually the same. Slight local vari- 
ations are due to Brownian motion of particles within 
their cages, but no rearrangements occur over the du- 
ration of the experiment. Due to the reproducibility of 
the experiment, it is reasonable to average the images 
"before" and the images "after" to reduce the variabil- 
ity caused by Brownian motion. The resulting difference 
image is shown in Fig. |31[d) and emphasizes the smoothly 
spatially varying displacements of the colloids. It can be 
seen that the applied force is not exactly on the x-axis of 
the image; this is dealt with in the analysis below. 

To quantify the images shown in Fig. we perform 
particle image velocimetry (PIV) on the pairs of "before" 
and "after" images. This method is frequently used in ex- 
perimental fluid mechanics, and does not depend on iden- 
tifying or tracking individual particles^. A small window 
in the first image is taken, and cross-correlated with the 
same size window in the second image. By moving the 
second window around in the second image, we find which 
piece of the second image is best correlated with the piece 
from the first image. The shift required for this maxi- 
mum correlation is a displacement vector reflecting how 
the particles have moved between the two images, and in 
particular represents the displacement vector for the cen- 
ter of the window. We use a window size that roughly en- 
compasses two particles, although our results are not sen- 
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FIG. 5. (a) Displacement field based on data shown in 
Fig-Htb). The arrows indicate displacements of the colloidal 
particles, (b) Residual displacement field after subtracting 
off the fit to Eqn. [2] The arrows are magnified by a factor 
of 5; in reality, the longest displacement vectors in panel (b) 
are 0.3 /im. The central region near the magnetic bead is 
removed for clarity. For both panels, the circles indicate the 
initial and final positions of the magnetic bead, which moved 
from right to left, and are drawn to scale. The scale bar is 10 
microns long. The data correspond to Figs. [S] and [B] = 0.49 
and Fniax = 0.29 nN. Note that a displacement vector is cal- 
culated for every pixel in the raw images; here only every 6th 
vector is drawn. 



sitive to this choice. The technique is merely correlating 
the images and the particles provide contrast to help this 
work. A typical displacement field is shown inFig.lSIa), 
corresponding to the images shown in Fig. |ll[a,b). For 
the PIV analysis, we use the individual raw images such 
as Fig. mja) which leads to Fig. Sfb), rather than the av- 
eraged images which lead to Fig. IDJd). After computing 
the PIV analysis for each individual pulse, we average 
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FIG. 6. Rescaled displacement vectors as a function of 6; 
compare with Eqn. [2] The points are the data and the solid 
line is the fit to the equation. The data correspond to Figs. |4] 
and[Sl^a), using only data with r > ro = sa/s. For this fit, a 
was constrained to be 1/2. 



the PIV fields over the pulse sequence for a given F^ax 
and (j) to do the subsequent analysis. 

The smoothly varying appearance of the displacement 
field seen in Fig. ^a.) suggests trying to fit the strain 
field to a simple functional form. As noted previously, 
the response of the magnetic bead to a constant force 
is a simple linear function, and so a natural choice is to 
treat the colloidal suspension as a homogeneous elastic 
medium. In such a medium, the strain field u around a 
point force F applied at the origin is given by^: 



1 1 



8nE 1 - cr 



{3-4a)F + h{n- F) 



(1) 



where E is the Young's modulus, a is the Poisson ratio, 
and h is a unit vector pointing away from the origin. 
Using n = X cos 9 + y sin 9 and F = Fx, the equation can 
be rewritten as: 



1 - cr 



[{7 -8a + cos 29)x + (sin 29)y] 

(2) 

which highlights the key spatial dependence of the strain 
field: it decays as 1/r, and has a periodic dependence on 
29, due to the symmetry of the problem about the x-axis. 
9^0 corresponds to the direction of the force. 

To test this, we rescale the displacements and Uy 
(measured from PIV) by r. This collapses the data rea- 
sonably well, as shown in Fig. |6l Here the data are plot- 
ted as a function of 9, showing the characteristic modu- 
lation in Eqn. [51 The solid line in Fig. [S] is a fit to the 
equation. The amplitude of both fit curves is constrained 
by the model to be the same, which is in slight disagree- 
ment with the raw data, where u^r is somewhat larger 
in amplitude than UyV. This is seen in most of our data 
sets. The curve for u^r is vertically offset, and it can 
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FIG. 7. Young's modulus _E as a function of volume frac- 
tion (f). Due to an inadequately defined applied force, E is 
overestimated although this affects all points equally (by a 
multiplicative factor) and does not change the shape of the 
curve; see the text for a discussion. The inset shows the same 
data plotted as a function of 0c — with 0c = 0.64. The lines 
in the main plot and the inset are the fit to the data using 
E = Eo{(t>c - 0)"'^ with Eo = 0.4 Pa and /3 = 1.84 ± 0.40. 
The symbol size indicates the uncertainty. 



be seen from Eqn. [2] that the magnitude of the offset is 
related to the Poisson ratio a. 

The fit has several parameters. First, the direction for 
= is chosen to be the average direction of all of the 
displacement vectors, to correct for the imperfect mag- 
net alignment. Second, the two physical parameters to 
the fit are a and E. A difficulty in determining E is that 
the true value of F is unknovifn: the model assumes a 
steady F whereas we apply a pulse. If the force was held 
at -Fmax, the magnetic bead would move with a veloc- 
ity and would not return to its original position, given 
the large values of F^ax we use^^. Fortunately, u scales 
reasonably well with -Fmax and thus leads to a consistent 
value for the Young's modulus E', even if its true magni- 
tude cannot be deduced from our fits. Third, we allow for 
the location of the origin (a; = 0, j/ = 0) to vary. It is not 
obvious if the origin should be at the starting position 
of the magnetic bead, the ending position, or elsewhere, 
especially given that the magnetic bead is a finite-sized 
disturbance and the model assumes a point-sized distur- 
bance. We adjust the origin so that the model has the 
best fit to the data; this typically puts the origin within 
0.3 /im of the starting position of the magnetic bead. A 
final parameter to our fitting algorithm is above what ra- 
dius To from the magnetic bead the fitting is conducted. 
Sufficiently close to the magnetic bead, its finite size be- 
gins to distort the strain field from the model. We fix 
''o — o,MB, and in practice our results are not sensitive 
to our choice. 

Despite the slight disagreements between the raw data 
and the model shown in Fig. \6\ overall the model is re- 
markably successful. For a few samples, we find that 
the Young's modulus E increases slightly with increasing 



^max, but more often we find E is independent of -Fmax 
and accordingly for each sample we average E over the 
different trials with different i^max- The resulting data 
are plotted in Fig. [3 and E increases by a factor of ~ 8 
as the glass transition is approached. Simulations and 
theory show that elastic moduli diverge near the jam- 
ming transition as B ~ (0c — with (pc ~ 0.64, the 
volume fraction of random close packing^^— . The ex- 
ponent /3 depends on the details of the interparticle in- 
teraction and which modulus is considered. The inset of 
Fig. [7] shows E plotted as a function of {(j>c — 4>) with 
behavior consistent with a power-law, although our data 
extend over only half a decade of {(jjc ~ </•)• Our exponent 
is /3 = 1.84 ± 0.40, similar to results for the bulk modu- 
lus of hard spheres (/? = 2) and shear modulus of hard 
spheres = 3/2)^^. 

Our values of E are quite similar to those found in 
a classic study of viscoelastic shear moduli of colloidal 
super-cooled liquids^^, although that is a coincidence. 
On the one hand, Eqn. [2] assumes the sample is in equi- 
librium for the applied force, which is certainly not the 
case. Using i^max overestimates E. On the other hand, 
our particles are 7.4 times larger than those of Ref J^, so 
our moduh should be smaller by a factor of 7.4'^ = 400. 
We can estimate the correct order of magnitude for our 
data from Fig. [51 using Eqn. [5] with our effective spring 
constant fc = 6.8 pN//im (setting this equal to u/F) and 
r — flMB = 2.25 fim. This gives us E ~ 0.72 Pa for 
4> = 0.55, suggesting that our data in Fig. [7] are overes- 
timated by a factor of O(IOO). Thus we are in plausible 
agreement with the data of Rcf.^^, in the high-frequency 
limit in particular which is most relevant for our quickly 
perturbed samples. 

An alternative comparison for E can be made with 
the theory of Schweizer and Saltzman, who developed 
an effective "free energy" for a hard sphere trapped in 
a cage^S. They construct the free energy F{r) as a 
function of the distance r from the cage center. For 
particle motion within the cage, they find an effective 
spring constant depending on as A: ~ fco exp(25.3(/)) 
with fco ~ 2.5 • 10~*kBT/a^ (where ks is Boltzmann's 
constant, T is the absolute temperature, and a is the 
colloidal particle radius). In our experiment, the mag- 
netic bead is larger than the surrounding particles by 
a factor of 1.45, so the effective spring constant experi- 
enced by our bead will be larger by 1.45^ = 2.1. Using 
4> = 0.55 and correcting for the bead size, their theory 
predicts k « 580kT/a^, as compared to our result of 
k = 6.8 pN/^m= 4000fcsT/a2. Our result is a factor of 
7 larger. Overall, given the approximations made by the 
theory and the uncertainties of the experiment, agree- 
ment within a factor of 7 is suggestive that the origin 
of the elasticity we observe is indeed the caging of the 
particles. 

The other key fit parameter in Eqn. [2] is the Pois- 
son ratio a. Over all values of i^max and (j) we find 
a = 0.50 ±0.08. Values of tr larger than 1/2 are unphysi- 
cal, so we conclude that our data show a = 1/2. Accord- 
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ingly, we fix this value and redo the fits to Eqn. [2l and 
the values of E that result are the ones shown in Fig. [7] 
and correspond to the fit curves shown in Fig. \6\ The 
physical meaning of ct = 1/2 is that volume is conserved 
during deformation: were this sample to be strained in 
one direction, the sides would contract sufficient to con- 
serve volume. This is plausible, as the sample itself is an 
incompressible fiuid with solid particles, and addition- 
ally one assumes the volume fraction stays homogeneous 
during simple deformations. 

As noted above, we allow the direction of the force to 
be a free parameter when performing the fit. This angle 
is fairly constant, with a standard deviation of only 4° 
between the different experiments. This variability likely 
reflects measurement error. 

The fit shown in Fig. [5] is not perfect, and some sys- 
tematic deviations from the fit can be seen. The differ- 
ence between the fit and the measurements is shown in 
Fig. \Elh). The displacement vectors are stretched by a 
factor of 5, and thus greatly exaggerate the difference. 
Nonetheless, this picture looks similar to the locally non- 
affine elastic behaviors seen in some simulation o^'^"" — 
and also images of "fioppy-modes," localized normal 
modes, and "soft spots" known to be present near 
jamming^"—. We stress that the majority of the total 
displacement field shown in Fig.[5l[a) is well-fit by Eqn. [21 

V. DECAY OF STRAIN 

A. Experimental observations 




0.0 0.5 1.0 1.5 2.0 

FIG. 8. Plots of the displacement of the magnetic bead as 
a function of time, after the force is removed, (a) Shows a 
linear-linear plot and (b) shows a log-linear plot. The values 
of Fmax are given in Table 1, with the largest initial displace- 
ment (red squares) corresponding to the largest force and the 
smallest initial displacement (purple pluses) corresponding to 
the smallest force. In (b), lines are fit to the initial data 
(t < 0.5 s) indicating decay time constants of 0.47 s, 0.38 s, 
0.37 s, 0.31 s, and 0.46 s (from largest -Fmax to smallest). 



After the force is removed, the magnetic bead moves 
back to its equilibrium position. Typical data of the mag- 
netic bead displacement as a function of time are shown 
in Fig.|8ja). Within our resolution, the magnetic bead is 
always in the initial position less than 10 s after it starts 
the return motion. 

Figure |8l[b) shows the data on a semilog plot, where 
straight lines would indicate exponential decay. While 
the initial portion of the data can be fit to straight lines, 
clear deviations are seen at longer times. The decay times 
found are 0.3 — 0.5 s but do not depend systematically on 
the initial displacement. Furthermore, some evidence of 
memory is seen. For example, the F^^x — 0.29 nN data 
(green triangles) go from x = 2.0 to 0.4 /im during the 
time interval t — 0.0 to 0.7 s. In contrast, the F^ax — 
0.75 nN data (red squares) go from x = 2.0 to 0.4 /im 
during the time interval t = 0.4 to 1.9 s, taking nearly 
twice as long to cover the same displacement. The noisy 
data seen in Fig. at small values of x are partly due 
to the uncertainty in determining x (±0.04 ^m). The 
a; = position is defined by an average at long times and 
so is more accurately defined. Within our resolution the 
positions shown in Fig. (H^b) have not quite decayed to 
X = over the time period shown. 

One trend is that samples with larger (p (closer to the 
glass transition) decay faster, as is shown in Fig.[Sl Given 



the nontrivial memory effects, it is not obvious whether 
to compare data at constant initial displacement or at 
constant Fmax- This distinction turns out to be unim- 
portant. Figureinija) compares two different volume frac- 
tions with the same initial displacement, and the data 
at larger (/> decay faster. Figure [9ljb) shows three differ- 
ent volume fractions with the same i^max, with the same 
trend, data for larger (p decay faster. This makes intuitive 
sense, as the elastic modulus is larger for larger (Fig. [7]). 
While viscous dissipation rises as the glass transition is 
approached^^, apparently the elastic contribution to the 
magnetic bead relaxation rises faster, resulting in a faster 
relaxation. 



While the position as a function of time does not ap- 
pear to decay exponentially (Fig. |8]) , plotting the data in 
Fig. [Tolas a function of -\/t suggests x ~ exp(— -^/t/to)- 
The value of to is slightly larger for larger initial dis- 
placements, although the data in Fig. [H are fairly par- 
allel within each panel, showing that Iq is not changing 
that dramatically. Iq is clearly larger for lower volume 
fractions cf). 
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FIG. 9. Relaxation curves for several experiments, demon- 
strating that the decay is faster for samples with higher (p. 
(a) Comparison of two samples with </!> as indicated, that have 
nearly the same initial displacement. For the (j> = 0.44 data, 
the force is -Fmax = 0.13 nN, and for the (f) = 0.49 data, the 
force = 0.29 nN. Two different instances are shown 

for the (j) = 0.49 data (triangles and pluses), (b) Comparison 
of three samples with the same force (i^inax = 0.29 nN) but 
different as indicated. 
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0.0 0.5 1.0 1.5 2.0 2.5 

FIG. 10. Displacement plotted as a function of y/i for (a) (f> = 
0.47 and (b) <j) = 0.49. The different symbols indicate different 
values of Fmax- The values of -Flnax are given in Table 1, with 
the largest initial displacement (red squares) corresponding to 
the largest force and the smallest initial displacement (purple 
pluses) corresponding to the smallest force. The straight lines 
indicate fits to ~ exp{~ y^tjto) . For (a), the values of to are 
0.50, 0.29, 0.31 s (top to bottom). For (b), the values of to 
are 0.23, 0.21, 0.19 s (top to middle). 



B. Model of relaxing bead 

To explain the stretched exponential decay process, we 
develop a model that treats the relaxation of stresses in 
the viscoelastic colloidal sample. Consider the relaxation 
dynamics of a magnetic bead of radius umb initially at 
a; = in a viscoelastic medium that is suddenly displaced 
at time t = by an amount xq due to an imposed force 
F. The bead was originally (at t < 0) in equilibrium 
and the displacement will result both in a force exerted 
on the external medium by the bead creating a stress 
field in the colloidal medium and a reaction force by this 
medium on the bead. The experiments show clearly the 
presence of both memory and a stretched exponential 
behavior for the bead relaxation. We show here that 
if the induced stress field relaxes in a diffusive manner 
then such behavior arises. The reaction force will thus 
in general be a function of both the applied force F that 
creates the inhomogeneous stress field around the bead, 
as well as time t due to diffusive relaxation of the stress 
field. This reaction force will tend to bring the bead 
back to its original equilibrium position due to the elastic 
forces exerted on the bead together with a viscoelastic 
drag force that will dissipate energy. Thus we can write 



the equation of motion 

mx = —m / ({t — s)x(s)ds — k{F,t)x, (3) 

where the first term on the RHS is the viscoelastic drag 
on the bead and the second term represents the elastic 
force on the bead in the presence of a relaxing force con- 
stant k{F,t). Because the motion is slow we can replace 
the viscoelastic drag by its viscous zero frequency limit. 
Namely defining 7 = ({t)dt we can rewrite Eqn. |3]as 

i ^ --fx - uj{F,tfx. (4) 

We estimate 7 — QiraMBri/rn using Stokes' law with 77 
being the effective viscosity of the colloidal medium, and 
we define a;(F, tf = k{F, t)/m. 

Our first challenge is to estimate k{F,t). We can see 
from the strain field u{r) induced in the colloidal medium 
due to the applied force F (see Eqn. [T]) that there exists 
a length scale ^0 over which the colloidal displacements 
will be greater than the typical colloidal particle radius 
a. This scale can be estimated as 

^o{F) « CF/iEa) (5) 

where E is the Young's modulus of the sample and the 
constant C « (l/167r)(H-a)(7-8cr)/(l-cr) « 0.2. Using 
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our largest F^ax = 0.75 nN and i? w 15 Pa, ~ 4a. 
Given that the medium is not perfectly elastic but rather 
viscoelastic, we argue that this region grows diffusively 
as the strain is dissipated in the surrounding medium, 
leadaing to a growing length scale £,{F,t) where 

aF,t)' ^^,{F)'+6Dt. (6) 

In Eqn. |6] -D is the diffusion constant characterizing mo- 
tion that can relax the strain, which docs not require 
cage rearrangements. We discuss D in more detail be- 
low. Using ^(F, t) we can estimate the typical strains 
in the colloidal medium induced by the bead. These 
strains e « x/^{F,t) reduce with time both due to the 
diffusive relaxation of the initial strain and the reduced 
imposed forces as the bead returns to its original equi- 
librium position. The associated elastic stresses in the 
colloid are then « Ex/£,{F,t). Thus we are now in a 
position to estimate the restoring force on the bead as 
k{F,t)x « 47rR'^Ex/^{F,t) or 

u:iF, tf « A^R^E/{mi{F, t)) = u^l/^il+t/r), (7) 

where w§ = Att E / {m^Q{F)) and r = £,o{Ff/QD. Sub- 
stituting Eqn. [7] into Eqn. U] then yields 

X = -7i - ijjlx/^J{l + t/r). (8) 

The relaxational behavior of the bead can now be 
found by solving Eqn. [5] subject to the initial conditions 
x{t = 0) = xo and x{t = 0) = 0. Though Eqn. [5] cannot 
be solved exactly, it has two limiting forms. For t <ti t 
Eqn. [8] reduces to it — ~jx — ujqX. This is the equation 
of motion for an linear oscillator with two overdamped 
modes 

x{t) = (xo/2){a+exp [~ja^t/2] 

+a_exp[-7a+i/2]}, (9) 

using a± = 1 ± ^1 — AlOq/j"^. In the limit t ^ t Eqn. [8] 
reduces to a; = ~jx — cjq jtx. In this limit we have a 
stretched exponential solution 

x(t) « (a;o/2)a+ exp (- ^ftjt^) (10) 

where to = 7V(^o^) = {'^^ I^^W^ Ikali^E''). Signifi- 
cantly, tg does not depend on the initial displacement x 
or the initial applied force F . 

We can compare these predictions to the experiment. 
As mentioned above, using Fmax — 0.75 nN and E « 
15 Pa, ^0 ~ 4a. The discussion in Sec. IIVI makes clear 
that neither Fmax nor this inferred E are the proper val- 
ues for Eqn. [SJ but on the other hand their ratio is what 
is needed to compute and it is precisely this ratio 
that is directly measured in the experiments of Sec. IIVI 
We estimate D as the short-time diffusion coefficient, 
D ss ksT /Q-nrja — 0.064 /^m^/s. This approximation 
using the dilute-limit value is imperfect due to hydro- 
dynamic interactions which reduce D at larger volume 



fractions^"—, but we are mainly seeking the right order 
of magnitude. Using this D and we find r « 100 s. 

The drag force acting on the magnetic bead is not due 
to the viscosity rj of the solvent (used to calculate D) but 
rather the effective viscosity of the medium, which is ~ 50 
times larger at these volume fractionsi^. To calculate to 
we use the more correct value of E estimated from the 
data of Fig. [2] as discussed in Sec. |TVl Using auB — 
2.25 /im and E — 0.72 Pa we get to = 1 ms. This is too 
small by a factor of ~ 200 from the experimental data 
(Fig. [To]). Likewise, given r « 100 s, we would expect to 
see the asymptotic (stretched exponential) behavior for 
\/(t) ^ 10 in Fig.[TOj that we see it at earlier time scales 
suggests that our estimate for r is too large. 

We thus reconsider the correct value of D. In our 
model, we assume D is the diffusion coefficient for strain. 
In practice, individual colloidal particles do not need to 
move significant distances for the strain to diffuse. Much 
as a dislocation can move rapidly through a crystalline 
lattice while individual particles stay close to their lattice 
sites, a slight motion of a particle (Ar < a) changes the 
strain over a neighborhood ^ a in scale. If we assume 
that particles diffusing a distance of a/20 is sufficient for 
the strain to diffuse a distance a, then D becomes 400 
times larger. This decreases r to 0.2 s and increases to 
to 0.4 s, bringing our model into more reasonable agree- 
ment with the data. The distance a/20 is smaller than 
the cage size (which is about a/3)i^. 



VI. CONCLUSIONS 

We have used magnetic beads to locally perturb a 
dense colloidal sample at volume fractions <j) < (pg, close 
to the colloidal glass transition. The magnetic beads have 
a linear relationship between the applied force and their 
displacement, and the strain field around the beads is 
well-described as that of a homogeneous elastic medium 
subject to a point force. The Poisson ratio is ti = 1/2, 
consistent with a sample that conserves its total volume 
when a stress is applied. Not surprisingly, the Young's 
modulus describing the elastic medium grows as the glass 
transition is approached. The growth is consistent with 
power-law in {(jjc — 4>), where (f)c = (f'rcp > 4'g- 

When the bead is moved away from its equilibrium 
position and the force is removed, we observed the sub- 
sequent relaxation to the equilibrium position. This 
relaxation behaves as a stretched exponential, x ^ 
exp(— (t/to)^/^). This agrees with a model that assumes 
the stress can diffuse away to infinity: thus, while the par- 
ticle is moving back to a; = 0, the effective spring constant 
acting on the particle is also diminishing. The experimen- 
tal time scales suggest that this diffusion is rapid, occur- 
ring faster than the particles themselves diffuse. This is 
likely due to the relatively small displacements of parti- 
cles needed to change the strain. 

We thank R. E. Courtland, S. A. Koehler, 
K. S. Schweizer, and M. Wyart for helpful discussions. 
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